Interplay between short and long-range entanglement in symmetry protected phases 
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We study a variety of questions related to entanglement in symmetry protected phases, especially those in- 
troduced in arXiv: 1106.4772 [1]. These phases are analogous to topological insulators in that they are short 
range entangled states with symmetry protected edge or surface states. We show that the now standard bulk- 
edge correspondence relating the entanglement spectrum to the gapless edge spectrum holds for these phases 
as well. We also consider the question of coupling these models to gauge fields or equivalently of introducing 
long range entanglement. We argue that this procedure yields models with perturbatively stable edge or surface 
states at a variety of interface types. The non-onsite nature of the edge symmetry plays an important role in our 
considerations. 



I. INTRODUCTION 

The main theme of this investigation is the nature of bound- 
ary states in various kinds of quantum matter. We now know 
of many examples of edge and surface states that are either 
completely robust (fractional quantum Hall states) [2-5] or 
are robust in the presence of certain symmetries (topological 
insulators) [6-9]. However, we still do not understand the 
general principles underlying the structure of edge states in 
2d and surface states in 3d, not to mention to complexity of 
purely boundary phenomena. 

The phases we study host a complex interplay of long and 
short range entanglement and promise to advance our under- 
standing of the physics of surface states arising in topolog- 
ical quantum matter. Such surface states are of experimen- 
tal interest given the edge modes in quantum Hall bars [4, 5] 
and the surface states of three dimensional topological insu- 
lators [8]. Many of the examples of fractionalized phases for 
which we have a current experimental candidate host gapless 
states, either in the bulk or at the boundary (see for example 
[4—8, 10, 11]). This is not surprising given the general visibil- 
ity of gapless states in experiment probes e.g. thermodynam- 
ics and transport. Of course, we are not at this time propos- 
ing that the phases we consider here are directly relevant for 
current experiments, indeed, our exactly solvable models are 
quite artificial. However, any advance in our understanding of 
gapless boundary states may be useful in future experiments 
given the prominence of these states in experimental examples 
of topological quantum matter. 

We have a variety of interests in this work. First, we wish 
to further understand the nature of symmetry protected edge 
states in systems with a unique ground state. Key questions 
here include the conditions for symmetry protection, classi- 
fication of edge states [9], and the role of perturbations e.g. 
disorder. Second, we wish to understand further the relation- 
ship between the spectra of physical edges and the entangle- 
ment spectrum [12]. Recall that the entanglement spectrum 
is essentially the spectrum of the reduced density matrix of 
a subregion in the bulk of the material. This connection is 
important not only because it provides a bulk method to iden- 
tify symmetry protected phases, but also because it suggests 
how the edge classification feeds back into the bulk classifi- 
cation. Third, we are interested in the extent to which short 
ranged symmetry protected phases can be strengthened by the 



addition of long range entanglement e.g. can we make these 
phases more robust in the presence of perturbations? We now 
elaborate briefly on this last point. 

In Refs. [1, 13] an exactly solvable model on a square lat- 
tice was given that has a unique ground state on any closed 
manifold and symmetry protected edge states on any open 
manifold. The gaplessness of the non-chiral edge states is pro- 
tected by an onsite Z 2 symmetry. Either the edge breaks the 
symmetry or the edge is gapless. While the symmetry is de- 
fined onsite i.e. the unitary generating the symmetry has the 
form U = Yl U r , the natural ground state is a multipartite 
entangled state associated with plaquettes of the square lat- 
tice. This tension leads to symmetry protection as codified by 
group cohomology and prevents the bulk Z 2 symmetry from 
being an onsite symmetry within the low energy edge degrees 
of freedom. Indeed, the global symmetry operation requires 
operators that act at each effective edge site as well as between 
each neighboring pair of edge sites. Importantly, this spatially 
extended symmetry structure is preserved under coarse grain- 
ing. However, this protection is only present as long as the 
onsite Z 2 symmetry is preserved. Enter long range topologi- 
cal order. 

Consider a long range entangled phase where the low en- 
ergy excitations are described by an emergent Z 2 gauge theory 
as in Refs. [14, 15]. Such a phase of matter need not possess 
any global symmetry, and we will assume no such symmetry. 
However, the gauge theory description implies an emergent 
Z 2 gauge symmetry in the deconfined phase. This is a some- 
what strange symmetry that cannot be broken by any local 
physical perturbation, hence it is a natural candidate to pro- 
tect the short ranged state described above. To clarify what is 
meant, consider a gaped system with global symmetry G. An 
isolated subregion R of that system can sometimes be labeled 
by a definite representation 7Z of the group G which is the 
total charge enclosed in R. This label is robust to G symmet- 
ric operations within R, but by adding G non-invariant per- 
turbations and by passing through a symmetry breaking state 
we can return to a symmetric state with any other charge VJ . 
Now consider an emergent gauge theory with gauge group G. 
We are sometimes still able to label isolated regions by a total 
charge (or more generally by anyon quantum numbers), but 
now no local perturbation can change this assignment. Phys- 
ically this is because there must exist an electric field labeled 
by 1Z that connects R with its environment. 
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We will argue that by coupling gauge fields to symmetry 
protected phases of matter we can produce systems with more 
robust edge and surface states. As we discuss in more detail 
below, such states may not be protected, that is guaranteed 
to be gapless by some symmetry (since there is no physical 
symmetry), but they can be robust in the sense that no small 
perturbation can destroy them. In this way we can trade sym- 
metry based protection for increased perturbative robustness. 

Independent of this work, gauging the symmetry of these 
symmetry protected phases has recently been considered in 
Ref. [16] as a device to prove results about edge states in the 
ungauged model. That work also considers interesting dual- 
ity transformations between symmetry protected phases and 
long range entangled phases, and we comment briefly on the 
additional insights into our problem provided by this duality 
technology. 

We will begin with some general motivational considera- 
tions before turning to a discussion of gapless boundary states. 
Then we will elaborate on the bulk-edge correspondence in 
various symmetry protected phases before turning to the ques- 
tion of gauging these phases. Our results include establish- 
ing a bulk-edge correspondence for the symmetry protected 
phases considered here and an analysis of the perturbative sta- 
bility of edge states in gauged models. We begin with some 
general motivating considerations in an effort to isolate the 
essence of symmetry protected edge states. 



II. MOTIVATING CONSIDERATIONS 

In many integer and fractional quantum Hall states, chi- 
rality and the bulk gap lead to an absolutely protected edge 
state, and in topological insulators we may have robust edge 
states so long as a certain symmetry is preserved. Examples 
of states protected by time reversal include the quantum spin 
Hall state in two dimensions [6, 7] and the topological in- 
sulator in three dimensions [8]. These states have likely been 
realized in HgTe quantum wells as reported in Ref. [17] and 
Bi2Te3 crystals as reported in Ref. [18] among others. We 
can also consider classes of states called fractional topologi- 
cal insulators as in Ref. [19-21]. Generally speaking, these 
states stand in relation to topological insulators in the same 
way that the integer and fractional quantum Hall effects are 
related. One feature of these states is that the gapless bound- 
ary modes may be stable in an RG sense to perturbations that 
break the symmetry as in Refs. [22, 23]. This added robust- 
ness is in keeping with the general intuition that fractional- 
ization enhances the perturbative stability of the state. Such 
intuition arises from the understanding that fractionalization 
goes hand in hand with long range entanglement, a property 
that is often robust to local perturbations. A general theme in 
the investigation of boundary states has been the topology of 
free fermion bandstructures, but this approach is not directly 
applicable to the symmetry protected states we consider. In- 
stead, as argued in Ref. [1], the low energy theory of the states 
in Ref. [1, 13] can be understood in terms of a 8 term associ- 
ated to a discrete symmetry. 

8 terms and associated Wess-Zumino-Witten terms for con- 



tinuous non-linear sigma models (NLSM) are well known in 
the theory of quantum magnetism (see Ref. [24] for a nice re- 
view). For example, consider a 2 + 1 dimensional bulk phase 
described by a gapped SU(2) non-linear sigma model with 
a 8 term. If the space on which the system is defined has 
a boundary, then the bulk 8 term, which is a total derivative 
when considering smooth field configurations, descends to a 
WZW term on the boundary of the manifold. The action of 
this system is 

S ~ / dtdx\x[d^Ud^}+ Swzw (2.1) 

with 

Swzw J dtdx du e^ vX 

xtr [(U^d„U)(U^d u U)(U^dxU)] . (2.2) 

The auxiliary space required to define the WZW term is here 
interpreted as the bulk 2 + 1 dimensional gapped bulk. Indeed, 
the integrand of the WZW term is nothing but the 2 + 1 di- 
mensional action for the 9 term action. The 1 + 1 dimensional 
SU (2) NLSM in the presence of a WZW term is an interest- 
ing system because, as was shown in Ref. [25], there exists 
a stable RG fixed point. Thus this system can enter a stable 
gapless phase at the boundary provided the SU(2) symmetry 
is not broken explicitly. 

For the purposes of our discussion here, there are two im- 
portant ingredients to this story. These ingredients arise from 
considering the two paradigmatic phases of a magnet, sym- 
metry broken and paramagnetic. The first ingredient is the 
fact that a continuous symmetry cannot break spontaneously 
in 1 + 1 dimensions. The second ingredient is the tendency 
of the WZW term to disfavor localized symmetric states e.g. 
a paramagnetic state. Since this second piece of intuition is 
rather more vague, we briefly elaborate on it. 

We can only write the WZW term as a local integral over 
spacetime provided we sacrifice the manifest symmetry of the 
Lagrangian. When we make a global symmetry transforma- 
tion on the Lagrangian, the manifestly symmetric terms are 
invariant without further thought, however, the WZW term 
changes by a total derivative. The action is still invariant pro- 
vided the space is closed, and this follows because we real- 
ize our WZW terms on the boundary of another bulk system 
(the boundary of a boundary is empty). The fact that the La- 
grangian is only invariant up to a total derivative is equiva- 
lent to the statement that the operator generating the symme- 
try transformation is not a product of terms acting on disjoint 
degrees of freedom. In other words, we cannot simply remove 
degrees of freedom from the boundary while keeping the sym- 
metry intact. In this sense, a localized symmetric phase of the 
NLSM is in tension with the WZW term, for example, the 
quantum state must not be smoothly connected to a product 
state if the symmetry is preserved. 

A simple example of this is phenomenon is provided by the 
spin- 1/2 Heisenberg chain as shown in Ref. [26]. Of course, 
the spin rotation symmetry is defined onsite, but the gapless- 
ness of this system may be understood in terms of an approxi- 
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mate "superspin" symmetry that combines the Neel TV and va- 
lence bond solid D order parameters into a single SU (2) ma- 
trix U — D+iN-a (we. impose the constraint D 2 + |iV| 2 = 1). 
The low energy theory is then a NLSM for the field U with a 
WZW term at level one. This theory possesses a stable RG 
fixed point and hence a gapless phase. Now in accord with 
our general comments above, it should not possible to realize 
the larger superspin symmetry in an onsite fashion. Define 
Mi = N t and M 4 = D and write the SU(2) WZW term as 

S WZW oc J dtdxdue IJKL e^ x 

x Mi «9 m Mj d u Mk d\ Ml ■ (2.3) 

This term may be cast as a two dimensional action at the 
cost of losing manifest symmetry. Set Mi = rrii and M4 = 
y/1 — m 2 ps 1 where the m,, should be regarded as infinitesi- 
mal and the symmetry acts by acts m, =— > m, + Si, Expand- 
ing the WZW term about this classical configuration gives 



Swzw oc / dtdxe^e ljk m i d tl mjd l/ m k (2.4) 



and so we see that 5 Swzw is 



SS 



WZW 



6iJ dtdxe^^d^m^rrik (2.5) 



which is only zero up to a total derivative. Microscopically the 
non-onsite nature of the symmetry follows because the VBS 
state consists of entangled neighboring spins while the Neel 
state is a product and no unitary formed from a product of 
onsite unitaries can connect these two states. We note that the 
appearance of an enlarged symmetry at the low energy fixed 
point, here an enhanced SU(2)l x SU(2)r symmetry, also 
appears to be a common feature of symmetry protected edge 
states. 

Now as we said, a discrete generalization of 8 and WZW 
terms were considered in Ref. [13] and argued to lead to a 
gapless edge. Unfortunately, we do not yet have much un- 
derstanding of this phenomenon from a field theory point of 
view. Nevertheless, early results in Refs. [1, 27, 28] sug- 
gest the following picture at least in d = 2: the edge theory 
is non-chiral but the symmetry is represented chirally on the 
edge modes. In other words, the left movers, say, carry the 
conserved charge while the right movers do not. Hence any 
backscattering terms that would gap the edge are forbidden by 
symmetry. Furthermore, the basic intuition described above 
about non-onsite boundary symmetry also carries over, and 
we will see below the non-onsite character of the boundary 
symmetry explicitly below. A possibility that must be consid- 
ered in the discrete case is spontaneous symmetry breaking at 
the edge (this is also an issue in continuous models in higher 
dimension). This still leads to "gapless" states at the edge but 
now associated with the order parameter direction in the sym- 
metry broken state. 



III. BOUNDARY STATES AND NON-ONSITE 
SYMMETRIES 

Let us now understand the role of the non-onsite action 
of the symmetry Q in protecting the low energy edge de- 
grees of freedom. Onsite means that Q may be written as a 
product Q = r Q (r) with (QM)t = (Q^y 1 for all r. 
Examples include rotation about the z-axis in spin systems, 
Q(0) = ® r exp (—i9S z ), and phase rotations in boson sys- 
tems, Q(9) = (^) r exp (—i8n r ) (n r is the boson number op- 
erator). In general, we would require that Q be defined as a 
sitewise product of unitary representations of some group. We 
will be interested in precisely those cases where the edge sym- 
metry is not onsite in the above sense (nor should it be short 
range equivalent to an onsite symmetry). 

Consider a bulk gapped symmetry protected phase with a 
physical boundary. Suppose that the ground state is invariant 
under the symmetry and fully gapped in the gap. Let us further 
imagine that at the edge there are a set of low lying states, not 
necessarily gapless, but lying much below the bulk gap. For 
the phases considered here, we can construct exactly solvable 
models where the edge is associated with a degenerate man- 
ifold of states. Adding perturbations yields the low energy 
manifold of interest. We wish to understand under what con- 
ditions this low energy manifold is guaranteed to form a gap- 
less state. Let Q represent the low energy part of the action of 
the bulk symmetry on the edge states. We will be much more 
explicit about this construction when we give some solvable 
models. The key issue will turn out to be whether Q has a low 
energy part that is short range equivalent, meaning equivalent 
under a short range unitary, to an onsite symmetry. 

To get some intuition, suppose Q is related via a finite depth 
circuit V to an onsite symmetry Q = V(&) reE Q r )V . 
Then clearly the edge can be in a gapped and short ranged 
state. Indeed, take any gapped onsite Hamiltonian He = 
'^2 ireE h r with [h r ,Q r ] = and consider its ground state 
\ip). We can always find such a Hamiltonian by slightly coarse 
graining the system (although translation invariance may need 
to be broken). The state V\ip) is clearly invariant Q, it is short 
range entangled, and is the ground state of a local Hamiltonian 
with a finite gap. What happens if we cannot relate Q to any 
onsite symmetry generator via a finite depth quantum circuit? 

Let us consider a simple Z2 example. We use the notation 
Z and X for the Pauli operators. Let Q be defined on N spin- 
1/2 degrees of freedom labeled by Z r = ±1 with action 



Q = Q]Jx r 



where 



Q = exp 



(1 - Z r )(l - Z r+1 ) 



(3.1) 



(3.2) 



This symmetry as defined is clearly not onsite, but it is also the 
case that it is not short range equivalent to an onsite symmetry 
(this will be clear later). As an additional comment, this sym- 
metry may at first seem quite special, but from a renormal- 
ization group perspective it is fairly general. In the extreme 
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infrared limit, we should be able to express any symmetry op- 
eration, even non-onsite symmetries, in terms of actions on 
neighboring spins. 

Let us now introduce a family of Hamiltonians H(g) that 
all commute with Q. To prove that no ground state of H(g) 
can be symmetric and short range entangled, we make two 
general assumptions: 

• [Invariance] The ground state \ip(g)) is unique, gapped, 
and invariant under Q. 

• [Local operators] There are local operators Z r satisfy- 
ing QZ r = —Z r Q (here this follows by definition of 
Q, but in general we will want to restrict attention to Q 
for which this is true). 

Using these assumptions we can prove the absence of sym- 
metric short ranged states. 

Let g* correspond to a state deep within a gapped short 
ranged entangled symmetric phase and write \ip(g*)) = 
J2z={z } c(z)\zi...zn). The invariance of \ip(g*)) implies the 
equation 

Q(z)c(-z) = c{z) (3.3) 

for all configurations z. An immediate corollary is that 
\c(z)\ = |c(— z)\. Since g* represents a point deep in a short 
ranged entangled state, we may write the wavefunction as 

c(z)=Hf(z r ,z r+1 ), (3.4) 

r 

and furthermore, given the extreme short ranged nature of the 
states, we need only consider a short chain of four sites (the 
minimal number). The extension to larger chains will be ob- 
vious. We also assume that f(z, z') = f(z', z). Now consider 
the configuration I1I2 — I3 — I4 where the subscripts indicate 
the site and we have periodic boundary conditions. Eq. 3.3 
implies that 

/(l,l)/(l > -l)/(-l,-l)/(l,-l) 

= -/(l, 1)/(1, -l)/(-l, -1)/(1, -1). (3.5) 

Clearly this can only be true if some of the /s are zero. 

Suppose /(l,— 1) = 0, then it follows that only config- 
urations of the form 1...1 and 0...0 have non-zero ampli- 
tude i.e. the state is a ferromagnetic cat state (in order to 
have invariance) or otherwise breaks the symmetry. Suppose 
/(1,1) = 0. The condition that \c(z)\ = |c(— z)\ immedi- 
ately implies that /(— 1,— 1) = and hence that the state 
is anti-ferromagnetic. The same result obtains if we suppose 
initially that /(— 1, —1) = 0. Thus we have shown that a sym- 
metric short ranged state is incompatible with the non-onsite 
symmetry Q. In fact, we know slightly more, since there are 
long range correlated states with c(z) given as above e.g. a 
state where c(z) = e~^" h '^/ 2 with hi the classical Ising 
Hamiltonian and j3 c the critical temperature. These kinds of 
states are also ruled out despite their power law correlations 
for some observables. 

Our proof is certainly more elementary than the sophis- 
ticated proof using matrix product states in Ref. [13], but 



we can also check that in other simple cases the same logic 
applies. We are currently developing this more elementary 
approach, which focuses directly on the non-onsite nature 
of edge symmetry, as an alternative to the group cohomol- 
ogy logic [29]. Note also that our wavefunction is equiva- 
lent to a matrix product state with matrices given by A™p = 

\/g(a, m)g(m, (3) (m the physical index and all indices run- 
ning over {1,-1}), and hence we make the same assumption 
for the form of the edge wavefunction as in Ref. [13]. One 
of the main virtues of our approach, which focuses directly on 
the non-onsite symmetry, is that it is generalizable to higher 
dimensions [29], a situation for which no formal proof of gap- 
lessness for non-onsite symmetries yet exists. 

Once we consider the possibility of spontaneous symme- 
try breaking, the generic statement is that the edge spectrum 
must always have low lying states. These low lying states 
may either be due to the presence of gapless edge modes or to 
the presence of degenerate ground states coming from spon- 
taneous symmetry breaking. Having recapitulated our under- 
standing of the protected low lying nature of the physical edge 
spectrum, we turn now to our second primary interest, the bulk 
entanglement spectrum. 



IV. ENTANGLEMENT SPECTRUM IN SYMMETRY 
PROTECTED PHASES 

Recall that the entanglement spectrum of a region R is de- 
fined using the reduced density pn = tr^(p) via the ansatz 
pit = e~ Ka . The spectrum of Kr is called the entanglement 
spectrum of R when the state p is the ground state. Note that 
the word "entanglement" is only appropriate when the state 
p is pure, otherwise not all the information encoded in the 
spectrum of Kb, describes entanglement. The entanglement 
entropy of R is simply the thermal entropy of the Hamilto- 
nian Kb at temperature one. More generally, an unentangled 
state p is described by an entanglement Hamiltonian with a 
single zero eigenvalue and an infinite gap to all other states 
i.e. e~ Kn is a projector of rank one. On the other hand, if Hb 
is gapless then the state p is highly entangled between R and 
R (Rs complement). Thus the entanglement spectrum pro- 
vides an interesting visceral measure of entanglement as well 
as permitting the computation of more precise measures such 
as the Renyi entropies S a (R) = (1 — a) -1 log (tx(p R )). Now 
we ask how the entanglement spectrum looks in a symmetry 
protected phase of the type we have been considering. 

Consider a bulk Hamiltonian Hb with gapped ground state 
|V>b). Suppose the Hamiltonian commutes with an onsite 
symmetry Qb- Again, onsite means that Qb may be writ- 
ten as a product Q B = (g) r Q (r) with (QW)t = (QM)- 1 
for all r. In general, we simply require that Q B be defined as 
a sitewise product of unitary representations of some group. 
Suppose also that the ground state is also invariant under the 
symmetry Q B - Q B \ip B ) = |V>b). 

If we now consider a subregion R, we find that the state of 
R is still invariant under the symmetry in the sense that 

PB = = ^FiiQEpQe 1 ) = QrPrQb 1 (4-1) 
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where Qr = (£) reR Q( r \ In other words, we have 
[Kr, Qr] = 0. Note that it is crucial for this result that 
we trace out degrees of freedom on which the symmetry acts 
in the standard way i.e. we don't split the symmetry ac- 
tion across the boundary. The boundary law for entangle- 
ment entropy plus the bulk gap suggests that Kr should be 
interpretable as a local Hamiltonian for a physical edge. Let 
|{m r }) be a basis for the states in R with finite entanglement 
energy (in practice we can add a high entanglement energy 
cutoff so that this basis is much smaller than the full basis for 
region R). Suppose that the labels {m r } are distinguished by 
local changes of configuration confined to the boundary of R 
as suggested by the finite gap and correlation length. We may 
thus ask how Qr acts on the states |{m r }). This action is 
given by a new matrix Q defined by 

Qii\m} = 'y y j Q mm '\m'). (4.2) 

Our key result is that Q mm i is non-onsite with respect to the 
local m degrees of freedom and is short range equivalent to the 
symmetry action on physical edge states. In other words, the 
entanglement eigenstates support the same non-onsite sym- 
metry action as the physical edge states. Since the entangle- 
ment Hamiltonian Kr commutes with a non-onsite symmetry 
and is local at the edge, our previous arguments imply that its 
spectrum is necessarily gapless e.g. the entanglement Hamil- 
tonian is truly gapless or spontaneously breaks the symmetry. 

We will show this non-onsite structure explicitly below for 
the Z2 model, however, the result is quite general and follows 
directly from the form of the fixed point wavefunctions in Ref. 
[1, 13]. Because those wavefunctions are direct products we 
can explicitly construct the entanglement boundary degrees 
of freedom m and the entanglement Hamiltonian is trivially 
local. The cohomology classification then implies that such a 
local Hamiltonian symmetric under Qr must have edge states. 
A similar result is immediately visible in the field theory for 
the non-linear sigma models we were considering above. The 
results of Ref. [30] imply that the entanglement spectrum of a 
relativistic NLSM, such as a NLSM with bulk 9 term, is equiv- 
alent in its universal properties to the physical edge spectrum. 
But as we already described above, the physical edge of such 
a gapped bulk phase with 9 term is a NLSM with a WZW term 
and hence (at least in 1 + 1 dimensions) is gapless. Thus the 
entanglement spectrum is also gapless like the physical edge 
spectrum and the action of the symmetry generator is again 
non-onsite. 

We can also use the recent results of Ref. [31] to give 
an alternate proof of the bulk-edge correspondence for some 
phases in the cohomology classification considered in Refs. 
[1, 13]. The proof proceeds exactly as in the previous para- 
graph. Ref. [31] provides relativistic field theories taking the 
form of Chern-Simons theories for several classes of symme- 
try protected phases. The results of Ref. [30] thus immedi- 
ately imply that the bulk-edge correspondence holds for these 
systems. We must only choose symmetric boundary condi- 
tions for the Rindler space cutoff in [30] to demonstrate the 
gapless edge modes. Symmetry breaking is also a possibil- 
ity, but as we described above the entanglement spectrum will 



still possess low lying states. Note also that this proof is quite 
different from the arguments above based on a lattice model 
formulation. Indeed, our arguments above apply to all the co- 
homology phases in Ref. [1], but the field theory approach 
still only works for a limited number of these models with 
known low energy descriptions. This is clearly a promising 
direction for further progress. 

We have seen that the entanglement Hamiltonian can nat- 
urally function as a local Hamiltonian for the system with a 
physical edge. Furthermore, the bulk gap in such an entangle- 
ment Hamiltonian goes to infinity away from the edge. This 
last requirement is crucial to obtain a boundary law for the en- 
tanglement entropy since any finite gap bulk system at finite 
temperature will give a volume contribution to the entropy. As 
long as we make this stipulation that the gap go to infinity in 
the bulk, it seems consistent to take the entanglement Hamil- 
tonian to basically be any local Hamiltonian for the topologi- 
cal state with edge. However, we must be quite careful here. 

If the actual Hamiltonian of the bulk material is transla- 
tion invariant with a translation invariant ground state, then 
despite the appearance of a non-translation invariant opera- 
tor Kr in the exponential, the state e" KR is actually trans- 
lation provided the temperature is kept at one. On the other 
hand, the generalized operator e~ l3Ka is not generically trans- 
lation invariant for f3 ^ 1. Nevertheless, these operators are 
interesting since they control the Renyi entropies S a (R) = 
(1 — a) -1 log (tr(p^,)). By sending j3 — > 00 we can project 
onto the entanglement Hamiltonian's ground state which we 
have argued should be the ground state of the physical sys- 
tem with edge. In particular, if the edge is gapless, like the 
edge of a quantum Hall state, then this entanglement ground 
state will possess long range correlations at the edge. This 
is another quite intriguing feature of the reduced density ma- 
trix in topological phases: where the original density matrix 
6 -Kr encoc jes a mixed state with decaying correlations, hid- 
den within this density matrix is a pure state with long range 
correlations. 

What happens when the entanglement Hamiltonian has a 
symmetry breaking ground state? Recall that even in a sym- 
metry broken phase, the spectrum is symmetric, but there is 
usually a breakdown of ergodicity that prevents the full man- 
ifold of symmetric states from being explored. However, the 
true ground state e.g. a cat state will also possess long range 
correlations. This in this case we see that the f3 — > 00 limit 
still produces a long range correlated state. We can use this 
observation to develop an interesting argument for a gapless 
surface in three dimensional systems. Recall that in two spa- 
tial dimensions a discrete symmetry breaking phase can sur- 
vive to finite temperature. If the surface of a three dimen- 
sional symmetry protected phase generically breaks the sym- 
metry, then the entanglement Hamiltonian must similarly have 
a symmetry broken ground state. Normally we think of a 
breakdown of ergodicity in such a state, but if we use the 
formal Gibbs state, then cluster decomposition is violated i.e. 
correlators do not factor at large distances. In particular, there 
are long ranged correlations in such a state, but long ranged 
correlations are incompatible with a gapped bulk state. Thus 
it must be that the thermal system described by cxp (—Kr) is 
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FIG. 1: The green dot represents a single site on the square lattice. 
The onsite Z-z symmetry acts by spin flip on each qubit and by a 
controlled Z operation on each neighboring pair of qubits. These 
two steps commute with each other and square to one. 

always in a paramagnetic phase. Either Kr has a symmetric 
low temperature phase (which is necessarily gapless) or the 
ground state of Kr breaks symmetry and the "entanglement 
temperature" must always be greater than the ordering temper- 
ature. If the entanglement Hamiltonian does have a symmetry 
breaking ground state, then we can detect this by taking the 
limit j3 — > oo for the operator e~° KR as above. Thus we see 
again that one simple characterization of the bulk of a sym- 
metry protected state is that the operator lim^oo e~ l3KR al- 
ways contains long range correlations due either to symmetry 
breaking or gapless states in the entanglement spectrum. 



TABLE I: Eigenstates of the Z2 symmetry Q 



states 


Q 


number 


name 


(1/V2)(|0000) ± |1111}) 


±1 


2 


l±> 


(1/V2)(|1000) ± |0111})(+trans.) 


±1 


8 


|1, ±)(+trans.) 


(1/V2)(|1100) ± |0011» 


Tl 


2 


|12, =F> 


(1/V2)(|0110)±|1001» 


Tl 


2 


|23,T) 


(1/V2)(|1010> ± |0101}) 


±1 


2 


|13,±) 




FIG. 2: The plaquette term of the CZX Hamiltonian involves 12 of 
the 16 qubits in a plaquette. The blue dashed lines indicates a a 
projector P involving those qubits. The red dashed box indicates the 
four qubits involved in the X4 operator. See the text for definitions. 



V. A 2D EXAMPLE: THE GAUGED CZX MODEL 

Having given a general picture of the protected edge modes 
at a physical boundary as well as the bulk-edge correspon- 
dence for the entanglement spectrum, we would now like to 
turn to the question of edge stability in the presence of long 
range entanglement. The simplest setting where we can carry 
out the process of coupling to gauge fields is in the non-trivial 
Z2 phase described in Ref. [13]. This phase is paramagnetic 
just like the large transverse field limit of the Ising model, 
but the phases are nevertheless distinct as demonstrated by the 
group cohomology arguments of Ref . [1, 13]. We will explain 
how to gauge this model, which is known as the CZX model 
(for reasons that will become clear), and then later discuss 
more general constructions. We will also give a preliminary 
discussion of the edge structure at the end of this section. 

Let us begin with a review the ungauged CZX model. The 
CZX model is defined on a square lattice with four qubits per 
site (see Ref. [13]). These qubits carry an onsite Z2 symmetry 
whose action is given 

Q = {X®X®X® X)(CZ ®CZ®CZ® CZ) (5.1) 

where X is the single qubit Pauli and CZ is a controlled Z 
gate acting on pairs of qubits. The symmetry structure is 
shown in Fig. 1. Note the peculiar "extended" structure of 



the symmetry action. In the computational basis the action of 
Q is 

Q\abcd) = (-l) ab+bc+cd+da \a+lb+lc+ld+l). (5.2) 

A simple computation shows that Q 2 = 1. It is convenient to 
diagonalize Q in order to classify onsite states as Z2 charged 
or uncharged. The eigenstates of Q are shown in Table 1 . 

The Hamiltonian of the CZX model is a sum of plaquette 
operators H p acting on twelve of the sixteen qubits involved 
in a plaquette: 

Hczx = J^-^p- (5-3) 

p 

The structure of the interaction is shown in Fig. 2. The 
operator X4 is 

X 4 = |0000) (1111| + |1111) (00001 (5.4) 

while the projecters P are 

P= |00)(00| + |11)(11|. (5.5) 

Focusing on a single four qubit site, we must consider opera- 
tors of the form 

Oabc.a'b'c' = \abc)(a'b'c'\ ® I (5.6) 
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FIG. 3: Here we show a physical edge in the CZX Hamiltonian 
together with the effective one dimensional description above it. 
The blue dashed lines represent qubits constrained to the subspace 
{|00}, |11}} due to a bulk projector P (from the plaquette term im- 
mediately below the edge). These two states are mapped to the ef- 
fective site represented by the red dot with states |0) = ] 00) and 
|1) = 1 11). As the figure indicates, in terms of these effective low 
energy edge states, the Z2 symmetry is not onsite. We must act with 
X on the tilde variables at each red site but also with CZ operators 
between pairs of sites. 

where the identity / acts on the outer corner. In fact, when 
dealing with the CZX Hamiltonian we need only consider op- 
erators with a = a', c = c', and b = b' + 1 (all arithmetic 
is mod 2). We label these operators O a bc which is short for 
O a bc,ab+ic- The action of Q on these operators is given by 

QO abc Q = (-l) a+c O a+lb+lc+1 . (5.7) 

A single plaquette term in the CZX Hamiltonian may be reex- 
pressed as 

tt _ _ ST r>W a k\ 

n V — 2-~i axba-^y a 2 ba 3 KJ a 3 ba 4 U a 4 bai \->-°) 

a\a2a>3CL&b 

where the (r) superscript labels the site. When acting with 
Q (1) <3 (2) <9 (3) <2 (4) on this operator we find that all the phases 
cancel and the terms are simply permuted about. This es- 
tablishes the symmetry of the CZX Hamiltonian under global 
symmetry transformations. We also note that the different pla- 
quette terms multiply to zero H p H p > = 

The ground state of the CZX Hamiltonian is beautifully 
simple. We simply place the four inside corner qubits of each 
plaquette, those in the red dashed box of Fig. 2, into the state 

|0000) + |1111) (5.9) 

and tensor over all plaquettes. Although the ground state is 
unique on a torus, the system has non-chiral edge states so 
long as the Z 2 symmetry is preserved. The extended nature 
of symmetry action on the edge, as mentioned in the introduc- 
tion, is shown in Fig. 3. The application of group cohomology 
ideas allows to prove the following result: the ground state of 
the edge cannot be written as a matrix product state with finite 
bond dimension in any phase where the symmetry action is 
preserved and described by a non-trivial element of the group 
cohomology. This leaves two alternatives: either the bond di- 
mension is infinite (gapless edge) or the symmetry is broken. 

Unfortunately, it seems that in the Z 2 non-trivial phase the 
edge must be fine tuned to be gapless. In other words, the 
edge generically breaks the symmetry. We describe this in 
more detail below. Note that this implies that two tunings 



are necessary to reach a gapless edge: we must impose the 
symmetry and then further tune one relevant operator to zero 
at the edge. This defect is not present in the Z^ models we 
briefly consider below or more generally. 

We now wish to add strings to the model so that the result- 
ing phase has long range topological order as in Refs. [14, 15]. 
We can always rewrite the operators O a bc in terms of the diag- 
onal states in Table 1 . This is convenient if we want to gauge 
the symmetry. To each link rr' of the square lattice we asso- 
ciate a Z 2 variable r* r ,. Consider a local Z 2 transformation 

0(r) = nW (,) ) lr ( 5 - 1() ) 

r 

specified by {x r } (x r € {0, 1}). Under this transformation we 
take r* r i — > r* r , (-l) Ir+v . This establishes the t z variables 
as gauge fields. t x is analogous to an electric field operator, 
and we interpret t x = ±1 as the absence or presence of a Z 2 
field line or "electric string". Using the variable r* r , we may 
"parallel transport" any quantity transforming under Q^ r > to 
a quantity transforming under Qv), More generally, using 
strings of the rs we can parallel transport a Z 2 charged quan- 
tity from one site to any other site. To construct the Hamilto- 
nian of the gauged model, we first modify the plaquette oper- 
ators of the CZX model. We express the plaquette operators 
in terms of the states in Table 1 for the four sites surround- 
ing a plaquette. This is a tedious but straightforward exercise. 
Then, for each term in the resulting expression, we append 
factors of Tf r , so that the resulting operator is only charged 
under the local symmetry of a single site in the plaquette, say 
the upper left corner. Because the original operator was in- 
variant under the global transformation, the new operator will 
be invariant under all local transformations. 

As an example of the first step in this procedure, consider 
the operator Oqoo on site r. This operator is 

Oooo = |000)(010|(|0)(0| + |1)(1|). (5.11) 

In terms of the states in Table 1 it may be written 

i(|3,+) + |3,-))((3,+| + (3,-|) 

+ i(|13,+) - |13,-))«13,+| - (13,-|). (5.12) 

The ± signs dictate how the states transform under so 
that 

Q {r) Q ooQ {r} = ^(|3,+) - |3,-))«3,+| - (3,-|) 
+i(|13,+) + |13,-))((13,+| + (13,-|). (5.13) 

Effectively, + states are Z 2 neutral while — states are Z 2 

(r) 

charged. To parallel transport O u00 so that it transforms under 
Q( r ' instead of Q( r > we modify it as follows 

i(|3,+)+r r V|3,-))((3,+| + (3,-|r r V) 
+^(|13,+) - r r V|13,-))((13,+| - (13,-|7£.<JL14) 
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This exercise must be carried out for the all O a bc operators. 

The full Hamiltonian is a sum of three terms. The first term 
is the modified plaquette term of the CZX model. The second 
term is a "magnetic" term for the gauge fields given by 

p rr'£p 

where the sum is over all plaquettes p and the product is over 
the four links in a given plaquette. The third term is an "elec- 
tric" constraint term that takes the form 

where the product is over all sites next to site r. This term 
simply says that the number of electric field lines entering 
site r must be equal to the Z2 charge of that site (mod 2). 
One may easily check that all three terms commute with each 
other. The only non-trivial checks involve the electric term, 
but commutativity follows because the other two terms either 
move a charge and create a string segment or create a closed 
string. 

The ground state of this model is simple enough. We take 
the Z2 plaquette state described above and tensor it with a 
state of the r variables given by t z = 1 on all links. We then 
project onto the constraint given by the electric term. This pro- 
jection attaches Z2 charges to the ends of open strings. The 
resulting state is a sum over all string states with strings re- 
quired to end at sites with non-trivial Z2 charge as dictated by 
the bosonic wavefunction. 

The model we have defined represents the extreme decon- 
fined limit of gapped Z2 charged bosons coupled to a Z2 
gauge field. To obtain the "true" gauge theory, we should im- 
pose the electric term as an actual constraint instead of ener- 
getically. We take this route below for convenience, but the 
we should bear in mind that this is not necessary since the 
gauge structure is totally robust. The low energy theory is 
then deconfined Z2 gauge theory coupled to a Z2 topological 
NLSM with 9 term as discussed above i.e. a gauged topologi- 
cal NLSM. At an interface between two deconfined phases, 
one with the 9 term and one without, a lower dimensional 
WZW term will appear. As an aside, in the high energy lit- 
erature there are subtleties associated with gauging a WZW 
term, but we have not yet understand the role of these obstruc- 
tions in our setting. We already noted that the edge theory 
tends to spontaneously break the Z2 symmetry. However, as 
we will argue in more detail below, the emergent gauge struc- 
ture prevents us from adding terms that explicitly break the Z2 
symmetry within the deconfined phase. Thus it appears that 
while a gapless edge still requires fine tuning, that tuning is 
less than in the ungauged model since the symmetry is now 
effectively unbreakable. Furthermore, the ungauged 2d mod- 
els with Zn (N > 2) symmetry (and potentially more general 
groups) have a stable gapless phase over a finite parameter 
regime. 

At low energies this model is essentially equivalent to the 
original gapped boson model plus gapped Z2 vortices. But 
importantly, because the local Z2 symmetry is an emergent 



gauge symmetry, no physical perturbation can break the Z2 
symmetry. Alternatively, stability of the phase implies the ex- 
istence of a renormalized gauge charge that always generates 
local symmetries. So long as such a renormalization proce- 
dure cannot change the emergent cohomological structure, the 
special edge states should remain. For example, the "global" 
gauge transformation given by x r = x leaves all gauge fields 
invariant. This transformation has the same spatially extended 
action at an edge as before, and since this structure is robust 
under renormalization, fluctuations of the gapped gauge fields 
should not be able to modify it. We now discuss a few differ- 
ent types of edges in these theories. 

Consider first a boundary between two deconfined phases, 
one with symmetry protected short range topological order 
and one without. Since nothing is happening to the gauge 
field, and since the symmetry cannot be broken by terms in 
the physical Hamiltonian, we expect the system will have ro- 
bust non-chiral edge states. However, we emphasize that these 
states are not protected. It is also possible for these edge states 
to spontaneously break the gauge symmetry and thus enter a 
Higgs phase. In this case, the gauge theory is locally screened, 
the emergent symmetry is lost, and the protected edge states 
may disappear. 

Now consider a boundary between the exotic deconfined 
phase and vacuum. Here the question of edge states is more 
delicate. We might expect them to become gapped due to 
proximity to a confined state since the emergent symmetry 
is effectively destroyed in such a phase. For example, some- 
times a confined state is smoothly connected to a Higgs phase 
of the gauge theory, and as we argued above, proximity to a 
Higgs condensate will destroy the edge states. However, this 
argument is heuristic only since the vacuum, or really, open 
boundary conditions in a lattice model, is not obviously the 
same as a confined state. Let us ask what happens in the lat- 
tice model, first in the ungauged model and then with strings. 

If the CZX model is given a boundary by removing all the 
plaquette terms above some line, the model now has a large 
ground state degeneracy. The projectors coming from the pla- 
quette term just below the boundary effectively set neighbor- 
ing boundary spins equal as shown in Fig. 3. The resulting 
degenerate manifold is already an indication of gaplessness, 
but to see a dispersion we must add extra perturbations to 
the edge. These perturbations must be consistent with the 
action of the edge symmetry shown in Fig. 3. The sim- 
plest such perturbation is an Ising coupling of the for Z r Z r i 
which commutes with the symmetry. Another possibility is 
X r + Z r -iX r Z r+ i. The unusual structure here is a conse- 
quence of the extended nature of the symmetry. For example, 
we already see that we cannot simply add a strong transverse 
field. We must also couple in an additional three spin term. 
This is precisely the unusual edge symmetry preventing us 
from reaching a short ranged symmetric phase. 

On general grounds we know that the edge Hamiltonian 
cannot be tuned to point where the ground state is a gapped 
symmetric state. This leaves two possibilities. Possibility 1: 
the edge possesses a gapless phase analogous to the contin- 
uous group case. Possibility 2: the edge always breaks the 
symmetry unless we fine tune. Possibility 1 requires the ab- 
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sence of all relevant operators in the gapless theory. In the 
conventional Ising model the transverse field can tune through 
a phase transition, however we have seen this term is not al- 
lowed on its own. Suppose the theory does possess a relevant 
operator. It must be the case that adding this relevant operator 
to the fixed point action leads to symmetry breaking regard- 
less of the sign of the coupling. For example, either sign of 
the ZZ term breaks the symmetry (either ferromagnetically or 
anti-ferromagnetically) in the absence of the transverse field. 
However, we do not expect that this behavior is generic, for 
example, the Z^ models we consider below do possess a gap- 
less phase at the edge. This is plausible given that a U(l) 
symmetry cannot break spontaneously and that Zm mimics 
U(l) for sufficiently large N. 

Suppose we similarly truncate the Hamiltonian of the 
gauged CZX model. Let us further suppose for simplicity 
that we impose the electric term as a hard constraint. This 
holds even at the edge where we must modify the nature of the 
electric constraint term since we are now missing the vertical 
gauge field. Everything still commutes and all electric strings 
remain tensionless. Let us first set t z = 1 on all links. This 
product state is the ground state before projection. Clearly the 
boson model has exactly the same spectrum as before. For ev- 
ery state of the boson sector in the gapless edge manifold we 
may project onto the constraint subspace to produce a physi- 
cal wavefunction. This attaches strings to the boson charges, 
but as we already said, these strings are tensionless. Hence 
the energy is not changed upon projection and all states re- 
main exactly degenerate. We note in passing that the local 
Z2 transformations look quite strange at the edge where the 
global symmetry is realized in a spatially extended fashion. 

Now consider perturbations to the edge. We can take sym- 
metric operators in the boson sector that will descend to sym- 
metric terms acting on the edge degrees of freedom e.g. of the 
forms we discussed above. We then gauge these operators as 
before to produce gauge invariant perturbations. Furthermore, 
these perturbations reduce to the usual symmetric forms dis- 
cussed above when we set t z = 1 on all links. Such terms will 
not gap the edge unless they cause it to spontaneously break 
the gauge symmetry. In the ungauged model we could add lo- 
cal terms, charged under the Z 2 symmetry, that would gap the 
edge states, but because this perturbation is charged we must 
attach a Wilson line to render it gauge invariant. Since this 
Wilson line must end on another gauge non-invariant operator, 
the gauge structure forces Z2 charged perturbations to come 
in pairs. Each such pair will be invariant under a global Z2 
transformation and will be local. If the Wilson line connects 
two perturbations on the boundary we obtain just another term 
of the type considered above. However, the Wilson line could 
also stretch into the bulk. These terms actually look danger- 
ous, as if charge is leaking off the edge, but this is an illusion 
since charge cannot escape into the gapped bulk. 

We have argued that in the non-fluctuating gauge field limit 
the gapless edge states necessarily remain. However, there 
is the question of what happens when the gauge fields fluc- 
tuate. Of course, these fluctuations are gapped so we might 
already guess that they are innocuous. Generally we think 
of a deconfined phase as one where we can smoothly set the 



gauge coupling to zero. However, what could such a fluctu- 
ating gauge field do? It will not confine any charges by as- 
sumption. It can certainly renormalize the edge Hamiltonian, 
but we have already discussed such perturbations. Since the 
string tension term hr x flip r z , we may analyze its effects 
for small h by studying the boson model in a slowly varying 
background field. Such a background field may not be trans- 
lation invariant, but according to Ref. [13] the edge states 
are robust to disorder so long as the symmetry is preserved. 
Simply put, we expect that sprinkling a few t z = — 1 bonds 
throughout the edge will not disrupt the gapless states since 
the symmetry is still preserved and the edge states may tunnel 
through short barriers. If we computed thermodynamic quan- 
tities, say, by summing over configurations in a path integral, 
then because each configuration has edge states, the gapless 
edge states will be visible in thermodynamics even with the 
gauge fluctuations. 

We have so far assumed, following Ref. [13], that static 
symmetry disorder will not localize the edge states. One piece 
of evidence in favor of this conclusion is the extended nature 
of the edge symmetry. No matter how we group sites at the 
edge, we cannot reduce the symmetry to an onsite form. In- 
deed, the proof of the absence of a matrix product state in Ref. 
[13] does not rely on translation invariance, so in this sense the 
edge cannot be in a short range entangled state unless the sym- 
metry breaks even in absence of translation symmetry. Indeed, 
our argument for the lack of short ranged states compatible 
with the symmetry can also be generalized to translation non- 
invariant situations. Nevertheless, the fate of the disordered 
model remains non-trivial, for example, the disordered spin- 
1/2 Heisenberg chain flows to an infinite randomness fixed 
point and while it remains highly entangled on average (hence 
forbidding a matrix product description), the structure of the 
excitation spectrum is quite unlike that of a conventional CFT 
or a z 7^ 1 scale invariant theory. Note also that for the spin- 
1/2 chain dimerization can immediately lead to a gapped sym- 
metric phase, albeit one that breaks translation invariance. Of 
course, it is now known that fermionic topological insulators 
have unusual localization properties as discussed, for exam- 
ple, in Ref. [32]. Clearly the effects of disorder are an inter- 
esting target for future studies. 



VI. GENERALIZATIONS IN VARIOUS DIMENSIONS 

We have given a detailed discussion of our ideas in the con- 
text of the CZX model, but they are more general. It is pos- 
sible to construct an exactly solvable boson model realizing 
every symmetry protected phase discussed in Ref. [1], and 
since the symmetries of these models are always defined on- 
site, there is no obstruction to gauging the symmetry. Thus 
in general we must consider G gauge theories coupled to G 
NLSMs with topological 8 terms. We hasten to add that this 
may not be the most general possibility. The models are of 
the Kitaev-Levin-Wen type with a non-fluctuating gauge field 
coupled to gapped charge matter and are exactly solvable. In- 
cluding gauge fluctuations is not expected to modify the situ- 
ation given the bulk. An interesting question arises as to what 
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exactly these phase are? Since they possess long range entan- 
glement and anyonic quasiparticles, and since there are a finite 
number of such phases with a given number of quasiparticles 
(see Ref. [14, 15, 33]), we expect that these phases are ac- 
tually of the type considered in Ref. [15], at least regarding 
their bulk properties. Indeed, the classification of Ref. [34] 
appears to coincide with the classification in Ref. [1]. This 
expectation is confirmed in Ref. [16] where it is shown that 
the bulk quasiparticles of the gauged CZX model correspond 
to those of the doubled semion model. 

For the purposes of the NLSM, we have argued that the 
main effect of the emergent gauge structure is to prevent ex- 
plicit symmetry breaking terms from appearing. Although we 
discussed the case of Z 2 , the gauging procedure and our ar- 
guments are completely general to any group G with a sym- 
metry protected short ranged phase. The effect of gauging the 
G symmetry is that explicit symmetry breaking operators are 
forbidden in the deconfined phase. The only remaining way to 
avoid gapless edge states is if the edge spontaneously breaks 
the symmetry. We know this behavior does not occur for con- 
tinuous groups, but these tend to confine in 2 + Id without 
a Chern-Simons term or gapless matter. For discrete groups, 
the Z2 case may require fine tuning, but more generally we 
expect a finite parameter regime in which the edge remains 
gapless. In other words, we expect the absence of symmetric 
relevant perturbations. Once we gauge the discrete symme- 
try explicit symmetry breaking perturbations are not allowed, 
thus the edge will be perturbatively stable to all deformations. 
We note again that the case of continuous symmetries is more 
complicated since these gauge theories often confine in two 
dimensions or are gapless in higher dimensions, furthermore, 
gauging a continuous symmetry on the edge can change the 
central charge of the edge theory potentially removing gap- 
less degrees of freedom, where as gauging a discrete symme- 
try (orbifolding) does not reduce the central charge [35]. 

We can also generalize our models to three dimensions. 
However, the exactly solvable models are extremely compli- 
cated, so we will not consider them explicitly. We begin by 
pointing out that the group cohomology classification is dif- 
ferent in three dimensions, for example, there is only one dis- 
tinct Z 2 short ranged phase. On the other hand, the groups 
Zn x Zn do have non-trivial short ranged phases. Thus we 
can consider an ideal Hamiltonian with surface states pro- 
tected by a Zjy x Z^ symmetry. We comment on the sur- 
face states and their robustness below. More generally, we 
can consider solvable models for any group G whose group 
cohomology indicates non-trivial symmetry protected states 
in three dimensions. We may also gauge the symmetry G 
using an identical procedure as in two dimensions. Given 
an ungauged Hamiltonian one can immediately write down 
a gauged Hamiltonian which is also exactly solvable. 

Regarding the ungauged edge model, we still have the 
NLSM story involving a boundary WZW term that is sup- 
posed to render the edge gapless. Unlike in the two dimen- 
sional case, there is no formal argument from, say, tensor net- 
work states that the surface state cannot be gapped, symmet- 
ric, and short range entangled. The intuition from Ref. [1] is 
simply that a gapped, symmetric, short range entangled phase 



cannot possess such a "non-local" term in its effective action. 
On a more practical level, the unusual surface symmetry ap- 
pears to prevent the addition of operators that could drive the 
system into a "paramagnetic" phase. Thus we see at least three 
possibilities: 

• the surface is gapped by spontaneous symmetry break- 
ing, 

• the surface is long range entangled because it is gapless 
(unless we fine tune the dispersion to zero), 

• or the surface is gapped and long range entangled be- 
cause it forms a topological state. 

Option three was not available for a one dimensional edge. 
There will be a gapless phase so long as all relevant operators 
are forbidden by the emergent symmetry and the WZW term. 
If there is no gapless phase, then again we find the intrigu- 
ing possibility that the surface either breaks the symmetry or 
enters a topologically ordered phase. 

Our direct approach to proving absence of symmetric short 
ranged states, as exemplified by our analysis of the CZX 
model above, can be generalized on a case by case basis, and 
we are currently working on a formal argument in three di- 
mensions and beyond. Although no formal proof along the 
lines of the matrix product state argument of Ref. [13] has 
been given for higher dimensional surface states, we believe 
this is likely a technical challenge. We believe the proper set- 
ting to formulate such an argument is provided by the multi- 
scale entanglement renormalization ansatz (MERA) (see Ref. 
[36]). MERA circuits of finite depth should be able to repre- 
sent any gapped, symmetric, short range entangled state. We 
are currently trying to better understand the action of the rele- 
vant non-onsite symmetry operators on MERA states in order 
to provide a proof of the conjecture in Ref. [1]. 

Let us conclude this section with a somewhat different ex- 
ample in three dimensions that may be divided into a "matter" 
sector and a "gauge" sector. Suppose that the gauge sector 
does not experience any qualitative change in its low energy 
structure over the whole space while the matter sector has a 
domain wall separating a symmetry protected phase from a 
trivial phase. We further assume that the symmetry protected 
phase has gapless boundary states at such a domain wall in the 
absence of the gauge field. In such a situation, we generically 
expect that the gapless states on the domain wall will remain 
intact even at finite gauge coupling. 

In fact, this situation is already realized experimentally, af- 
ter all the physical electron is charged under the electromag- 
netic gauge field, and although we often treat electromag- 
netism as a background field, it does fluctuate in real sys- 
tems. Such a realization is provided by the interface between 
Bi2Te3 in its topological phase and vacuum where the gauge 
field is the physical electromagnetic field. Now experimen- 
tally these surface states are seen in ARPES (see Ref. [18]), 
so we must conclude that within the resolution of ARPES the 
gapless states remain even at finite gauge coupling. 

But what is the gauge field is not free outside the material? 
Returning to the case of fractionalization, a convenient ex- 
ample is provided by the 3 + 1 dimensional topological Mott 
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insulator proposed in Ref. [37] and further analyzed in Ref. 
[38]. This state may be understood by decomposing the elec- 
tron operator as c a = bf a where & is a boson that carries the 
electron charge and f a is a fermion carrying the electron spin. 
This decomposition has a 17(1) redundancy given by b — > e l6 b 
and f a — > e~ l6 f a and hence b and / are coupled by a U(l) 
gauge field. Let us consider four distinct phases of this system 
based on whether the boson is in a superfluid or Mott insulat- 
ing phase and whether the fermion is in a trivial or topological 
band insulator. If the boson is condensed and the fermion is 
trivial then we obtain a trivial insulator of electrons (I), and 
if the boson is condensed and the fermion is non-trivial then 
we obtain a topological insulator for electrons (TI). In both 
these phases the U(l) gauge field is screened by the boson 
condensate and effectively disappears from the spectrum. If 
the boson is in its Mott insulating phase, then the gauge field 
is no longer screened and we have two kinds of fractionalized 
phases depending on whether the /s form a topological (FTI) 
or trivial insulator (FI). 

We now describe the six possible interfaces between these 
states. We assume throughout that time reversal and the elec- 
tromagnetic U(l) are unbroken. We know that at an I/TI inter- 
face there are gapless electronic states. We also know, based 
on our discussion above, that at an FI/FTI interface there will 
be gapless states of the / electrons. At an I/FI interface we do 
not expect any surface states since the boson is simply con- 
densing and screening the gauge field. Similarly, at a TI/FTI 
interface we also do not expect to find surface states since the 
bandstructure of the / fermion is not changing. At a FTI/I 
interface we expect gapless electronic states (rather than gap- 
less states of the / electrons) since the condensed boson in the 
/ phase will, via the proximity effect, screen the / fermions 
with bosons to produce electrons near the interface. Finally, at 
an FI/TI we also expect electronic surface states. We should 
point out that those interfaces where more than one thing is 
happening e.g. the boson is condensing and the fermion is 
changing its bandstructure are more sensitive to microscopic 
details since we can always arrange for the changes of phase to 
occur in microscopically different regions. For example, the 
proximity effect depends exponentially on distance and hence 
our arguments can depend on the precise length scales of the 
interface. A clean way to handle this is to approximate the 
interfaces in question with double interfaces where at each 
interface only one component (boson or fermion) changes it 
state e.g. I/FTI w I/FI/FTI or I/TI/FTI. 

One very important interface type we have not yet consid- 
ered is the interface with vacuum. In many respects the vac- 
uum behaves like a trivial insulator, but there are subtle dif- 
ferences. Because the vacuum contains no electron charge, 
the boson b should certainly not be regarded as condensed in 
the vacuum. On the other hand, the emergent U(l) gauge 
structure is certainly absent in the vacuum since it emerges 
from the collective dynamics of electrons and there are none 
present. We may refer to the vacuum as a confined phase of 
the gauge theory in the sense that physical electrons injected 
into the vacuum will not fractionalize into b and / particles. 
So the question arises as to what surface states, if any, exist 
at an interface with vacuum. We can make a few statements 



right away. A vac/I interface has no protected states while a 
vac/TI interface does. It is also believed that a vac/FI interface 
has no protected edge states, but what about a vac/FTI inter- 
face? Suppose we model this interface as a vac/FI/FTI inter- 
face, that is, an interface where, as we emerge from the bulk, 
the fermions change their state before the bosons and fermions 
go to zero density and the gauge theory confines. Then we will 
have protected surface states of the / fermions coming from 
the FI/FTI interface; furthermore, nothing obviously singular 
happens as the width of the FI region shrinks to zero. On the 
other hand, if we model the vac/FTI interface as a vac/TI/FTI 
interface then we have surface states of electrons. In the first 
case the fermion first changes its bandstructure while in the 
second case the boson first condenses, but the second case is 
special because as we already noted, the boson is certainly not 
condensed in the vacuum. Hence something singular is hap- 
pening as the width of the TI region shrinks to zero, namely 
the condensate (b) which is converting the surface state of / 
fermions into a surface state of c electrons is potentially van- 
ishing. We can certainly have a stable surface phase where b 
is not condensed. If the condensate vanishes, then we might 
again expect to recover surface states of the / fermions. This 
is indeed what was found in Ref. [38] where the boson is sim- 
ply uncondensed and the gauge field is shown not to destroy 
the gapless surface states (although it does modify the details 
of the surface). 

Another perspective on the question of surface states and 
boundary conditions is provided by the bag model of quark 
confinement introduced in Ref. [39] (see Ref. [40] for an easy 
introduction). This model, originally developed in the con- 
text of confinement in QCD, supposes that quarks and gluons 
are confined into regions called "bags". These bags may be 
roughly understood as deconfined bubbles within a larger con- 
fined vacuum. The spectrum of quarks and gluons confined to 
such a bag with appropriate boundary conditions then repre- 
sents an approximation to the spectrum of QCD. Furthermore, 
the bag is not static but may fluctuate in spacetime leading to 
a variety of interesting dynamical effects. Besides modeling 
hadrons, the bag model has also been applied to situations 
like quark matter deep inside neutron stars where a decon- 
fined quark phase gives way to a confinded phase of hadrons 
as the surface of the neutron star is approached. Below we will 
analyze our U(l) toy model in the context of the bag model 
and show that within this framework there are indeed surface 
states at a vac/FTI interface. 

Let us work in 3 + 1 dimensions and consider a bag con- 
taining a Dirac fermion ip and a £7(1) gauge field a M . We will 
assume that the bag surface is fixed e.g. at the edge of our 
material and does not fluctuate. The gauge field couples to a 
conserved current given by j M = ip^ip. Our conventions for 
the 7 matrices are 



with ct p = (1, a 1 ) and = (1, — a 1 ). The requirement that 
no charge leave the bag is given by the constraint n^j^ = 
where n is a unit normal to the bag boundary. This constraint 
may be written explicitly as %pij ■ nip — 0. Using the fact that 
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(7 • n) + = —7 • n for spacelike n and (17 • n) 2 = 1, we have 
the freedom to choose 17 • nip = ±ip at the inner boundary of 
the bag. This also implies that ijjip = at the bag boundary. 

Now suppose that in addition to the bag boundary, the ip 
fermion mass also changes sign at the bag boundary. With- 
out the bag there would be a surface at this interface, but is 
this zero mode consistent with the bag model boundary con- 
ditions? Outcome. Let x 3 be the direction normal to the local 
bag boundary (we have in mind a very large bag). The Dirac 
equation with the varying mass m(x 3 ) = m sgn(x 3 ) is 

ij^d^ip = m{x 3 )ip (6.2) 

and we look for a solution with d^ip = for fj, ^ 3. However, 
we see immediately that a solution can be found if ip satisfies 

ij 3 ip = sgn(a; 3 )-!/> where ij 3 = ij ■ n. Note that if ij 3 tp = 
±ip then we have iptp = -0+7°^ = 'ip ± ij 3 ip, but this last 
expression satisfies 

± ipij 3 ip = ±tp + (-ij 3 )j°i{,> = zp(i7 3 -(/))+7 ^ = -ifaP 

(6.3) 

which implies that ipip = in such a situation. Thus the sur- 
face mode is compatible with the bag boundary conditions. 

If we now include the gauge field, we must specify further 
boundary conditions. The original bag model boundary con- 
ditions of Ref. [39] are n^F^,, = which ensure that the 
bag is always a charge neutral object. One heuristic argument 
for these boundary conditions comes from the dual supercon- 
ducting model of confinement. Just as magnetic monopoles 
are linearly confined in a conventional superconductor (by 
vortex lines), the electric charges in a dual superconductor 
are similarly linearly confined. In this model the dual vari- 
ables F^ v = ^e A1 „Ao-^ lA<T satisfy conventional superconduct- 
ing boundary conditions e lu ,\ n v F*' a i.e. that the tangential 
dual electric field and normal dual magnetic field vanish at 
the boundary, but these conditions imply n^F^. In a gauge 
where A 3 — Owe find that the condition F^i = may be sim- 
ply reduced to d^Ai = at the boundary. Furthermore, the 
condition F30 = reduces to 83 A = 0. The resulting model 
of surface modes coupled to a U(l) gauge field with these 
boundary conditions has been analyzed in Ref. [38] where it 
was found the surface state survives (those authors choose a 
slightly different gauge, but we do not expect this to seriously 
affect the results). 

We have seen that topological Mott insulators can have sur- 
face states at an interface with vacuum (provided time reversal 
is unbroken). On the other hand, these states are not abso- 
lutely stable, that is they can be gapped without closing the 
bulk gap. One possibility is that b condenses near the surface 
in which case there will remain edge states but of electrons 
and not / fermions. Another possibility is that the / fermions, 
which sit in a finite density metallic-like state on the sur- 
face, could spontaneously pair and condense. In this case the 
surface behaves like a superconducting layer for the internal 
U(l) gauge field. We see that there are numerous possibilities 
and the complete analysis is complex. If the electromagnetic 
U(l) and time reversal are unbroken then we certainly have 
surface states of some type, but in general the surface states 
break up into different perturbatively stable classes. Which 



surface states are realized in a given model then depends on 
the microscopic details of the interface. 

VII. DISCUSSION 

Here we list our main results and conclusions. 

• [Non-onsite symmetries] The non-onsite nature of the 
symmetry prevents, in a rather direct way, the existence 
of short ranged symmetric ground states. 

• [Bulk-edge correspondence] The entanglement spec- 
trum of a large bulk region reproduces the universal fea- 
tures of the physical edge spectrum e.g. it is gapless or 
spontaneously breaks symmetry. 

• [Fractionalized phases] For any short range phase in 2d 
or 3d protected by an onsite symmetry G, we can con- 
struct a stable fractionalized phase with an emergent G 
gauge symmetry so long as the G gauge theory has a de- 
confined phase. We can further construct exactly solv- 
able models for these phases. 

• [Robust edge states] In some cases these fractionalized 
phases possess perturbatively robust gapless states at 
a boundary between two distinct deconfined phases as 
well as at a vacuum boundary. 

The last point deserves further comment. We can distin- 
guish between edge states that are totally robust, edge states 
that are totally robust in the presence of a symmetry, and those 
that are perturbatively robust. The latter category refers to 
edge states that can be removed by sufficiently strong pertur- 
bations but which are described by an RG fixed point with 
no perturbatively relevant directions i.e. for small deforma- 
tions. The exemplar of the first category are the quantum Hall 
edges while the second category describe systems like topo- 
logical insulators. Our conjecture apparently can only refer to 
the third category since it seems that a weakly coupled Higgs- 
type phase at the boundary would always remove any putative 
boundary states. 

We have we argued that by weakly gauging the symmetry in 
a symmetry protected phase, we can trade symmetry protec- 
tion for perturbative stability. The symmetry protected edge 
cannot be gapped without breaking the symmetry (protection) 
yet there are relevant perturbations that produce a gap e.g. 
those that break the symmetry. The gauged symmetry pro- 
tected edge can be gapped (no physical symmetry protects it) 
but it can be perturbatively stable e.g. the gauged Zjs; (N > 2) 
case discussed above. Since such gauged edges are not pro- 
tected, the question naturally arises as to what extent they can 
be regarded as stand alone theories in a lower dimension with- 
out the bulk. This questions would be interesting to study in 
the future, so here we only offer a few brief comments. 

From the point of view discussed at the end of section 2, 
namely that the edge may be understood as carrying the sym- 
metry chirally, it seems unlikely that such state could exist on 
its own. For example, consider the U(l) states discussed in 
Refs. [27, 28]. These states all have an even Hall response for 
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a background gauge field coupled to the conserved U(l), but 
this implies that the edge current suffers an anomaly e.g. it is 
not conserved in general. The explicit process which removes 
charge from the edge corresponds to flux insertion in the bulk. 
The existence of the Hall response for the background £7(1) 
gauge field implies that letting the £7(1) gauge field fluctuate 
leads to a gapped topological state in the bulk (confinement 
of pure £7(1) gauge theory is avoided). However, certainly 
the gauged chiral edge theory is inconsistent without the bulk. 
Note that the bulk theory is somewhat similar to the chiral spin 
state discussed in Refs. [41, 42]. In particular, in Ref. [42] it 
is argued that the degrees of freedom associated with the £7(1) 
edge modes are removed by the gauging process, and since we 
have no additional SU (2) symmetry it appears that the chiral 
edge states originally carrying the £7(1) charge are completely 
removed from the spectrum. This leaves the intriguing possi- 
bility that the remaining non-charged modes with the opposite 
chirality could survive on the edge. This possibility deserves 
further study. 

The phases we have constructed are, at least as far as the 
bulk is concerned, identical to previously known long range 
entangled phases as argued in Ref. [16]. This connection is 
important because the models constructed in Ref. [15] are not 



generally expected to have protected edge states, so this corre- 
spondence agrees with our expectation that the gauged theo- 
ries we consider do not have protected boundary states. Nev- 
ertheless, perturbatively stable gapless boundary modes are 
quite interesting in their own right. We also reiterate the point 
that while the bulk may be identical, there appears to be sub- 
tle questions related to boundary conditions especially in the 
duality arguments. It is possible that while there is only one 
bulk doubled semion phase, there are multiple renormaliza- 
tion group stable boundary conditions, some of which support 
gapless edge states. General considerations along these lines, 
although for gapped boundaries, have recently been given in 
Ref. [43]. We should also point out that some of the mod- 
els we have constructed are amenable to study via quantum 
Monte Carlo and hence the question of edge states should be 
answerable numerically for those models. 
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